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We report a neutron diffraction and muon spin relaxation /iSR study of static and dynamical 
magnetic properties of BaCo2V2 0s, a quasi-one-dimensional spin-chain system. A proposed model 
for the antiferromagnetic structure includes: a propagation vector /caf = (0,0, 1), independent of 
external magnetic fields for fields below a critical value H C (T). The ordered moments, of 2.18 {is 
per Co ion, are aligned along the crystallographic c-axis. Within the screw chains, along the c 
axis, the moments are arranged antiferromagnetically. In the basal planes the spins are arranged 
ferromagnetically (forming zig-zags paths) along one of the axis and antiferromagnetically along 
the other. The temperature dependence of the sub-lattice magnetization is consistent with the 
expectations of the 3D Ising model. A similar behavior is observed for the internal static fields 
at different muon stopping sites. Muon time spectra measured at weak longitudinal fields and 
temperatures much higher than Tjv can be well described using a single muon site with an exponential 
muon spin relaxation that gradually changes into an stretched exponential on approaching Tjst. The 
temperature-induced changes of the relaxation suggest that the Co fluctuations dramatically slow 
down and the system becomes less homogeneous as it approaches the antiferromagnetic state. 

PACS numbers: 75.10.Pq, 75.25.-j, 75.40. Cx, 75.40. Gb, 75.50. Ee, 76.75.+i 



INTRODUCTION 

One-dimensional (ID) spin chain systems with nearest- 
neighbor antiferromagnetic (AF) exchange interactions 
have attracted much attention in the last decades be- 
cause they display a rich variety of ground states with 
non-classical magnetic phenomena. [TH3] It is well known 
that an ideal ID AF spin system does not undergo long- 
range order (LRO) at finite temperatures due to strong 
quantum spin fluctuations. [4] But slight perturbations 
of the system by, for instance, weak interchain interac- 
tions or doping with non-magnetic impurities, can make 
the quantum critical state unstable and favors three- 
dimensional (3D) AF LRO. [5] 

In 2005 BaCo2V20g, a new member of the quasi-lD 
AF spin-chain family, was synthesized by Zhangzhen He 
and co-workers. [7H9] They succeeded in growing large sin- 
gle crystals of high quality. BaCo2V20s is an insulator 
with a very anisotropic magnetic susceptibility the 
easy axis of magnetization is along the crystallographic c 
axis. The temperature dependence of the magnetic sus- 
ceptibility x(T) displays a broad maximum near 30 K 
signaling the development of ID short-range order. [7] 

Crystal-field effects and spin-orbit coupling partially 
lift the degeneracy of the single ion Co 2+ (3d 7 ) multiplet 



resulting in a doublet ground state. Thus the spin-system 
may be considered as a collection of effective S — 1/2 
spin chains. [10 With lowering temperature, a rapid drop 
of x a t Tn = 5.4 K signals the onset of a long-range or- 
dered AF state at zero field. [7] Quantum fluctuations of 
the spins strongly affect the ordered state, as evidenced 
by the phenomena observed in external fields applied 
parallel to the chains, which results in an unusual field- 
temperature phase diagram (see Fig. 1). 

The magnetization curves M{H) are strongly non- 
linear. [TOj-Qj] External magnetic fields H above a critical 
value H C (T) suppress the 3D AF Neel state. At T = 
K, fioH c — 3.9 T. The field-induced suppression of the 
Neel state occurs in an unusual manner. Above 1.8 K, 
H induces a reentrant phenomenon from the long-range 
3D, Neel state into a ID quantum liquid state (order- 
disorder transition) . [T5] This transition was theoretically 
predicted first for gapped spin systems doped with non- 
magnetic impurities. [14] 

At lower temperatures, however, external fields H > 
H C (T) result in the development of an incommensurably 
modulated AF state. The existence of this novel phase 
was demonstrated experimentally by Kimura and co- 
workers. 15-18] They found for the modulated structure 
a propagation vector Um = (0, 0,1 + 5) with S tunable by 
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FIG. 1: Schematic representation of the Magnetic Field - 
Temperature phase diagram of BaCo2V2 08 for fields applied 
along the c-axis. Two magnetically ordered phases are dis- 
tinguished: (i) The Neel state, a 3D AF state at low fields, 
and (ii) an incommensurably modulated phase at higher fields 
(from Kimura and co-workers, Ref. 17). 



the external magnetic field. 

Although many interesting properties have been inves- 
tigated in BaCo2V20g, there are no reports on the details 
of the structure of the Neel ordered state, i.e., in external 
magnetic fields H < H C (T). In addition, there is very lit- 
tle information about the spin dynamics in this system, 
except for two NMR studies. [HI [20] We present here a 
detailed model for the magnetic structure, obtained from 
neutron diffraction measurements in powder samples and 
in single crystals of this material. We also report the re- 
sults of /iSR measurements below and above T?j. The 
latter shine light into the dynamics of the spin system. 

The physical properties of BaC^A^Os can be under- 
stood in the context of Ising spin chains with weak in- 
terchain interactions as described by the S = 1/2, XXZ 
model where a single chain in the presence of an external 
magnetic field H is described by the Hamiltonian 

H = jJ2(SiS!+x + <SfS* +1 + SVS! +1 ))- 9t i B s?h 

i i 

(1) 

From magnetization and ESR data it has been deter- 
mined for BaCo 2 V20s that the AF exchange interac- 
tion is J/ks = 65 K, the anisotropy parameter e 0.5, 
and g ps 6.2. The interchain interaction J' is of the or- 
der of J'/ J « 0.06. [11] In this model Neel order occurs 
at zero field with a non- vanishing gap. [5] The applica- 
tion of external magnetic fields, H > H c , restores the 
Tomonaga-Luttinger spin-liquid state with longitudinal 
and transversal incommensurate spin fluctuations that 



characterize ID Ising systems. [2"l |2"T] 

BaCo2V208 crystallizes adopting a tetragonal struc- 
ture with space-group symmetry /41/acd, room- 
temperature lattice parameters of a = 12.444A and 
c = 8.415A and 8 formula units per unit cell. All the 
magnetic Co ions are equivalent. They lie in the centers 
of CoC>6 octahedra, and have strong AF exchange in- 
teractions with Co ions in the neighboring edge-sharing 
CoOe octahedra. The Co 2+ spins, here referred simply 
as "spins", form a collection of screw-type chain struc- 
tures along the c axis. The screw chains are separated 
by nonmagnetic VO4 (V 5+ ) and BaC>4 (Ba 2+ ) tetrahe- 
dra, resulting in a quasi-one-dimensional structural ar- 
rangement of the spins with a strong Ising-type magnetic 
anisotropy. A representation of the crystal structure is 
given in Fig. 2. 
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FIG. 2: Representation of the basal plane projection of the 
crystal structure of BaCo2V208- This figure emphasizes large 
blue circles, representing Co ions with the rods joining them 
representing the n.n. Co-Co bonds. The small orange circles 
represent the O ions, the green and violet circles represent the 
Ba and V ions, respectively (colors on line). 



EXPERIMENTAL PROCEDURE 

Poly crystalline samples of BaCo2V20g were synthe- 
sized by a standard solid-state reaction method using 
high purity reagents of BaC0 3 (4N), CoC 2 04-2H 2 (3N) 
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and V2O5 (4N) as starting materials. Single crystals of 
BaCo2V2 0s (of dimensions of the order of 10 x 8 x 3 
mm 3 ) were grown by a spontaneous nucleation method 
from the same starting materials. Details of the mate- 
rial synthesis have been previously described. [71 [HI [H] 
Our material was characterized by X-ray diffraction and 
magnetic susceptibility. The crystal structure and the 
magnetic properties confirmed the conclusions based on 
previous data. [13] 

The magnetic structure of BaCo2V2 0s in zero exter- 
nal magnetic field was investigated by measuring neutron 
powder diffraction patterns at temperatures down to 1.8 
K using incident neutrons of wavelength of 2.4503 A in 
the Cold Neutron Powder Diffractometer DMC at SINQ, 
Paul Scherrer Institute (PSI), Switzerland. The sample 
consisted of 8 g of powdered BaCo2V2 0s material en- 
closed in a thin-walled V cylindrical container. On a sin- 
gle crystalline sample, we measured the temperature and 
field dependencies of some selected magnetic Bragg re- 
flections using the Triple Axis Spectrometer Rita-II and 
the Thermal Single Crystal Diffractometer TriCS both 
at SINQ, PSI. The latter with incident neutrons of a 
wavelength A = 1.18A. External magnetic fields, of up 
to fi H — 4.5 T, were applied along c axis. /iSR time 
spectra were measured on a polycrystalline sample, in 
the temperature range between 1.5 and 100 K in zero 
and weak longitudinal magnetic fields, using the Gen- 
eral Purpose Surface-Muon Instrument GPS at the Swiss 
Muon-Source of PSI. 
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FIG. 3: Bragg scattering profiles for BaCo2V208 at T = 8 
K (a) and T — 1.3 K (b). A pattern corresponding to the 
difference of intensities [7(1.3 K)— 7(8 K)] and revealing only 
the magnetic Bragg reflections is shown in (c). The solid lines 
represent the refinement of the data using a model described 
in the text. 



RESULTS AND DISCUSSION 
Magnetic structure 

In Figs. 3(a) and 3(b) we display neutron powder 
diffraction patterns above (8 K) and below (1.3 K) the 
Neel temperature, respectively. The magnetic Bragg re- 
flections, i.e., the difference of intensities [7(1-3 K)— 7(8 
K)] is displayed in Fig. 3(c). A refinement using the 
"FullProf" software package [55] yielded a very good fit 
of all the magnetic and nonmagnetic reflections. From 
the nuclear reflections we obtained for the lattice param- 
eters of the tetragonal cell at 8 K: a = 12.389 A and c = 
8.375 A. All the magnetic reflections, which appear be- 
low T/v, could be indexed in the tetragonal system with 
a wave vector Icaf — (0,0, 1). The best fit to the data 
are represented by solid lines in Figs. 2(a)-2(c), where 
the conventional Rietveld factors are 7? p = 6.08, 7? wp = 
6.13 and 7? oxp = 1.13. 

The refined spin structure is represented in Fig. 4, 
where the Co ions are depicted as circles and their bond- 
ings as the heavy solid lines joining them. Different co- 
ordinates along the c axis are indicated by the different 
fillings of the circles. This figure emphasizes the parts 
of the screw chains contained within the unit cell. The 



chains run along the c-axis, the ordered Co moments are 
aligned along the same direction and "u" and "d" rep- 
resent up and down spin directions. The spins are ar- 
ranged antiferromagnetically along the screw chains, but 
at a given basal plane they are arranged antiferromag- 
netically (ferromagnetically) forming zig-zag paths along 
the a (b) axis. The Bragg peaks are very narrow, reso- 
lution limited, thus, the 3D long-range character of the 
AF order is confirmed. From the refinement, the magni- 
tude of Co moments is estimated to be 2.18 /xb which for 
an effective spin 5=1/2 implies an effective g-factor of 
<7«4.4. 

Our neutron powder-diffraction data are not compat- 
ible with the value g = 6.2, obtained from high-field 
ESR and magnetization results. [7 Assuming the ESR 
5-value, the calculated intensities of the magnetic reflec- 
tions would be too large, by a factor of 2, to fit the exper- 
imental data. Therefore, we conclude that the Neel state 
involves reduced Co moments. This reduction is proba- 
bly due to quantum fluctuations expected in the 3D XXZ 
model. But the discrepancy of the (/-factors extracted 
from neutron scattering, high-field ESR and magnetiza- 
tion results have to be considered carefully since these 
techniques probe the spin system under different condi- 
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FIG. 5: The temperature dependence of the square root of the 
integrated intensity of the (3, 2, 0) magnetic reflection. The 
inset shows the magnetic field dependence of Bragg intensity 
of (2, 1, 0), (4, 1, 0), (3, 2, 0) magnetic reflections and the 
nuclear (2, 0, 0) as a function of magnetic field along the c 



FIG. 4: Representation of the magnetic structure of 
BaCo2V208 in the AF state. Co ions are depicted as circles 
and their bondings as solid lines. The quasi-lD screw-type 
spin chains run along the c-axis and the "up" and "down" 
spin directions are represented by "u" and "d" , respectively. 



tions. 

We have measured the temperature and field depen- 
dences of selected magnetic Bragg reflections on a single 
crystalline sample using Rita-II. A summary of these re- 
sults is presented in Fig. 5. The circles in the main 
panel represent the square root of the integrated mag- 
netic Bragg intensity yTe of the (3, 2, 0) reflection as 
a function of temperature. This quantity is proportional 
to the sublattice magnetization. Its functional form, i.e., 
the \Z7b(T) curve, follows the expectations for a 3D Ising 
antiferromagnet, here represented by the solid line. In a 
limited temperature range near Tjy (0.8 < T/T/v < 1), 
we find 
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(2) 



The fits to the data using Eq. (2) yielded a Neel temper- 
ature Tn — 5.4 K, as expected, and a critical exponent 
/3 = 0.28. This value is close to the predictions of the 3D 
Ising model, which has been estimated to be (3 w 5/16, 
with the exact value depending slightly on the approxi- 
mations used. For the 2D case: /3 = 1/8. 011] 

The inset of Fig. 5 shows the normalized integrated 
intensity of three magnetic Bragg reflections (2, 1, 0), (4, 
1, 0) and (3, 2, 0) as a function of magnetic field along the 
c axis. For comparison, we also plot the intensity of the 
nuclear Bragg reflection (2, 0, 0). The data indicate that 



the magnetic ordering is suppressed by external magnetic 
fields of 3.75 T along the c axis, at 2.4 K. This is in 
agreement with the results obtained from measurements 
of the specific heat in external magnetic fields. [8] [TT] We 
found no measurable field dependence for the magnetic 
propagation vector Uaf m the Neel ordered phase. 

In an effort to add more light to the nature of the phase 
transition into the spin liquid phase, we have searched for 
quasi-elastic scattering in the vicinity of the above men- 
tioned magnetic Bragg reflections. The searches were 
performed at temperatures near Tjv and zero external 
magnetic field using the thermal single crystal diffrac- 
tometer TriCS with a 2D detector, at SINQ, PST No 
clear quasi-elastic signals were detected. 



Muon spin relaxation and rotation 

To investigate the spin dynamics of BaCo2V20s in 
the paramagnetic state, we have performed /iSR experi- 
ments in zero-field ZF and in weak longitudinal external 
magnetic-fields LF.[57] Four typical /LtSR spectra above 
Tjy, measured in a field of 20 Oe are shown in Fig. 6. This 
weak external field effectively suppresses the contribution 
to the muon relaxation due to quasi-static nuclear dipo- 
lar fields. This leaves basically a dynamical relaxation, 
which for the case of a single relaxation rate is usually 
well described by a simple exponential function. [28 More 
complex dynamic cases are described using stretched ex- 
ponential functions. (29] [30] 

The LF-/xSR time spectra data above Tjv were fitted 
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FIG. 6: /J.SK time spectra measured in a weak longitudinal 
external magnetic field, of 20 G, in the paramagnetic state at 
four temperatures well above TV. The solid lines represent 
fits to the data using effectively a single component (see text). 
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FIG. 7: Temperature dependence of muon relaxation rate Alf 
(a) and stretched exponent n (b) obtained in a external lon- 
gitudinal field of 20 G in the paramagnetic state. The inset 
in (a) depicts the distribution of relaxation rates at T w 8K, 
where n — 0.5 (see text). 



using 

G z (t) = Aexp[-(A LF i) n ], (3) 

where A is the initial asymmetry and Alf is a parameter 
also denoted here as the muon longitudinal spin relax- 
ation rate. The best fits are represented by the solid lines 
in Fig. 6. In the paramagnetic state the relaxation (of 
dynamic character) is associated with the fluctuations of 
the Co moments. An stretched exponent n ~ 1 reveals a 
sharp distribution of rates, which effectively may be rep- 
resented by a single spin relaxation rate Alf for the entire 
spin system, while smaller values, n < 1, suggests a dis- 
tribution of rates which rapidly broadens with decreasing 
n. Only two parameters are needed to fully describe the 
distribution of rates: The (stretched) exponent n char- 
acterizes the width of the distribution and 1/Xlf is the 
time needed for the initial asymmetry A to decay to 1 [e 
of its initial value (independent of n). 

The results of the temperature dependence of Alf and 
n are summarized in Figs. 7(a) and 7(b). At high tem- 
peratures, T > 30K, the data are well fitted with an 
stretched exponent n « 1, signaling that the system is 
rather homogeneous. Below 30 K, with decreasing tem- 
perature, Alf (I 1 ) displays a substantial increase suggest- 
ing a rapid slowing down of the dynamics of the Co spin 
fluctuations. These are considered as in the "fast fluc- 
tuation regime" at all temperatures of these particular 
measurements, T > 8 K. At the same time the system 
becomes more inhomogeneous, as judged by the rapid 
decrease of n. At 100 K, n 0.9 but just above Tn, 
n ps 0.5. The increase of Alf with decreasing T is as- 
sociated with the increase of AF correlations among Co 
moments within the individual screw chains. This claim 
obtains support from previous findings based on results 
for x(T). A broad maximum in x(^) signals the onset 
of intra-chain magnetic correlations. [5] An illustration of 
the inhomogeneity of the system at low temperatures is 
given by the distribution of rates depicted in the inset of 
Fig. 7(a). 

Stretched exponential relaxation appears in a wide va- 
riety of phenomena, including spin glasses and AF 
systems with frustrated interactions. In such systems, 
the relaxation at high temperatures follows an exponen- 
tial form and, as the temperature decreases (towards the 
freezing temperature in spin glasses) , the functional form 
changes into a stretched exponential function[30 

G z (t) = Aexp[-(\t)^}. (4) 

Here the stretched exponent depends on the effective di- 
mensionality d. Qualitatively this general picture applies 
to our system with a reduced dimensionality, since al- 
ready at 8.2 K, Fig. 7(b) shows that n < 0.5, clearly 
below the expectations for a 3D case (n < 3/5). There 
is a clear trend of decreasing n with decreasing temper- 
ature. At low temperatures the tendency seems to be 



6 




O 1 1 1 1 1 1 1 

0.1 0.2 0.3 

Time (us) 



FIG. 8: fiSK time spectrum measured in zero field in the AF 
state at 1.6 K. The solid line represents a fit to the data using 
5 different components. 

n — > 1/3, i.e., towards the d — 1 case. But this trend 
changes abruptly upon ordering, where the effective di- 
mensionality of the system suddenly transforms into 3D. 

In Fig. 8 we present the ZF-/1SR time spectrum at 
T = 1.6 K, deep in the AF state. Here a spontaneous 
Larmor precession is revealed by the time-evolution of 
the asymmetry. The prominent structures revealed in 
this plot signal the presence of several incquivalent muon 
sites with non-vanishing internal static fields. This is 
characteristic of a magnetically ordered state. 

The muon spectrum below Tjv is well described by as- 
suming four components with different relaxation rates 
and different Larmor frequencies: 

4 

G z (t) = Aexp(-Ai) + V" Ai exp(-A»t) cos(27r^t + <f>), 
i=i 

(5) 

where the first term, a signal with vanishing Larmor fre- 
quency, is associated with muons whose spins are oriented 
along the internal static fields. This corresponds to 1/3 
of the total signal for our polycrystalline sample, i.e., 
A fa 1/3(A + A 1 +A 2 + A 3 + At). Therefore, the contri- 
bution to the asymmetry from muon sites with vanishing 
static magnetic dipolar field must be small, otherwise the 
fraction of A relative to the total asymmetry would be 
larger. The signals with non-vanishing Larmor frequen- 
cies: v\ = 5.7, V2 = 11.1, V3 = 18.4 and V4 = 23.9 
MHz, at 1.6 K, arise from muon stopped at sites with 
internal static fields of 0.42, 0.82, 1.36 and 1.76 kOe, re- 
spectively. The largest asymmetry is 0.12 corresponding 
to and is substantially larger than the other signals 
(A4 ^S> A\ ~ A2 ~ A3). Muon motion in the crystal 
may be ruled out by the lack of appreciable temperature- 
induced changes in the asymmetries A\, A 2 , A 3 and A4 
below Tn. 

The model described by Eq. (5) is plausible. For 
instance, the largest relaxation rate corresponds to the 
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FIG. 9: Temperature dependence of non-zero muon relaxation 
frequencies inferred from the zero-field data in the AF state. 
All the frequencies vanish at Tjv = 5.4 K. The solid lines 
represent fits to the data using the 3D Ising model. 

largest internal fields; A4 ~ 28 MHz. The other three 
have substantially smaller relaxation rates, of the order 
of Ai ~ A2 ~ A3 ~ 7 MHz. The lowest rate corresponds 
to the signal with vanishing Larmor frequency A < 0.8 
MHz. Although the model yields good fits of the data it 
has limitations; namely, it is difficult to identify the muon 
sites. We found several interstitial sites, (x, y, z) with 
x/a, y/b, z/c = {0, 0.5}, where the electric field vanishes. 
However, at these sites the calculated dipolar magnetic 
field also vanishes and, therefore, they are unlikely to con- 
tribute to the observed asymmetry. Muon sites where the 
dipolar magnetic field correspond to the inferred internal 
static fields are, for instance, near the centers of sides and 
faces of some BaC>4 tetrahedra. However, at these places 
the electric field does not vanish. We conclude that most 
likely the observed signals must arise from muon bound 
to oxygen ions. 

In Fig. 9 we display the temperature dependence 
of the non-zero muon Larmor frequencies inferred from 
the zero-field data. All the frequencies vanish as T ap- 
proaches Tjv- The solid lines represent fits to the data 
using the critical behavior v = vo(l — (T/Tjy) 13 ). Experi- 
mentally one finds f3 ~ 0.28, similar to the value obtained 
in the analysis of the neutron diffraction data, and close 
to the expectations of the 3D Ising model. Therefore, 
below T/v the data consistently suggest that the ordered 
moments of BaCo2V2 0s may be considered as a 3D Ising 
system. 

In summary, as BaCo2V20s approaches the antifer- 
romagnetic state with decreasing temperatures: (i) the 
Co moment fluctuations undergo a drastic slowing down 
and (ii) the system, as probed by the muon relaxation, 
becomes gradually inhomogeneous. These features re- 
flect the development of ID intra-chain AF correlations 
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and show that just above T)v the dynamics of the spin 
system may be thought of as a disordered Ising system 
with reduced dimensionality (probably ID). The phase 
transition at T/v may be thought of as the order-disorder 
transition of ID objects. Upon order the dimensional- 
ity of system changes into 3. In the AF state, near T/v, 
the temperature dependence of the internal static fields 
follows the expectations of a 3D Ising model with an ef- 
fective S = 1/2. The magnetic structure is characterized 
by a wave- vector Uaf = (0, 0, 1), independent of temper- 
ature and external magnetic field. The moments of 2.18 
Hb per Co ion are arranged AF within the screw chains, 
along the c axis, and in the basal planes they form fer- 
romagnetic zig-zag paths along one of the axis and an- 
tifcrromagnetic paths along the orthogonal (equivalent) 
direction. 
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